NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



NASA Technical Memorandum 79281 
Ciast-Ta-79281) Mianic ttsfoist of oaiaaiB 

UOLmiBB PIBBI COnPOSXTBS CIISI) 17 p 

nc I02/BP aoi csct iib 


63/28 


DYNAMIC RESPONSE OF DAMAGED 
ANGLEPLIED FIBER COMPOSITES 


C. C. Chamis, J. H. Sinclair 
and R. F. Lark 
Lewis Res^rch Center 
Cleveland, Ohio 


Prepared for the 

Winter Annual Meeting of the 

America.i Society of Mechanica I Engineers 

New "York, New York, December 2-7, 1979 


Mo-iiias 

Onelas 

86PS1 



1|Kl»^,¥V.-'h 




omMiic RsspomE or OMiA<an> amgleplibo pibbr composites 

C. C. CbaaiSf J. R. Sinclair ana S. P. Lark 

MASA-bcwia Resaarcto Canter 
Cleveland. Oitio 44135 


! 

I 


ABSTRACT 

An investigation was conducted to determine the effects of low 
level damge induced by nonotonic load, cyclic load and/or resid- 
ual stresses on tte vibration frequencies and daaqping factors of 
fiber composite angleplied laminates. Two different coaq>osite 
systems were studied - low modulus fiber and ultra high modulus 
fiber composites. The results obtained showed that the fre- 
quencies and damping factors of angleplied laminates made from low 
modulus fiber cooq;>osites are sensitive to low level damage while 
those made fr<» ultra high modulus composites are not. Also vi- 
bration tests may not be sufficiently sensitive to assess concen- 
trated local damage in angleplied laminates. Ai'.d furthermore, 
dynamic response determined from low-velocity im^ct coupled with 

the Fast Fourier Transform and packaged in a minicom^ter can be a I 

convenient procedure for assessing low-level damage in fiber com- t 

posite angleplied laminates. 

INTRODUCTION \ 

j 

A major concern in the fiber composites community has been the | 

static and dynamic response of slightly damaged (not visually dis- 
cernable) angleplied fiber composites. Coupled with tnis concern 
is the need to use a convenient mechanical test to identify both 

the effects and the extent of damage. In response to this concern , 

and need, an investigation was conducted at Lewis Research Center 

(LeRC) with the primary objective to assess the dynanic response 1 

of slightly damaged angleplied laminates. A secondary objective 

of the investigation was to identify a convenient dynamic test 

that can be used to assess this damage. The investigation was 

both experimental and theoretical in nature. 

J in the experimental part, angleplied laminates K45, 9 O 2 , 

(quasi-isotropic) were made from two composite systems, 
system was E-glass fiber/epoxy matrix (E-G/E) and the other 
was ultra high modulus graphite fiber/epoxy matrix (T-75/E) . ‘ 

These combinations of laminate configuration and composite system 
were selected to bound the magnitudes of elongation to fracture 




and the M9nitudea of laaination residual stresses. Speciaens 
froai these laainates were subjected to cyclic aechanical load 
tests until a change in tie stress-strain curve was cbservcd. 
Speciaens were also loaded aonotonically to fracture. After these 
tests, segaents froa the speciaens were i.'.bjected to low level 
iapact tests. Their dynaaic response (frequencies and daaping 
factors) were deterained using the Fast Fourier Transfora Method 
available in a ainicoaputer. The speciaens tested were ultrason- 
ically C-scanned prior to and after the cyclic testing. 

In the theoretical portion of this study, laainate analy**is 
was used to deteraine first ply failures including both laaination 
residual stresses and those due to aechanical loao. NASTRAM was 
used to deteraine the vibration frequencies cf the speciaens sia- 
ulating the test conditions. The aechanical properties required 
by NASTRAN for both undanaged and daaaged speciaens were generated 
using laainate theory. 


EXPERIMENTAL 

The experiaental part of the investigation consisted of fabri- 
cating, instruaenting, and testing the speciaens. 

Laainate Fabrication 

The laainates were aade at the Lewis Research Center using 
coaaercial Thornel-75 graph ite/PR288 (T-7S/E) and E-glass/lP03 
(E-G/E) prepregs. Individual plies of prepreg were stacked at the 
required fiber orientations in aetal aolds to fora the 12 ply 
angleplied laainates. A theraocouple was inserted into one end of 
each uncured laainate. For the T-75/E laainate the cold aold was 
placed in a hydraulic presa the platens of which had been heated 
to 450 K (350° F) . A contact pressure of O.lO MPa (15 psi) was 
applied on the aold. This pressure was maintained for 3 minutes 
after the laainate reached 311 K (100° F) ; then the pressure was 
increased gradually over a 2 minute period to 2.1 MPa (300 psi). 
This pressure was maintained for 2 hours after which the acid con- 
taining the cured composite was taken from the press. The lamin- 
ate was then removed and allowed to cool to room temperature.. For 
the E-G/E laainate, the press platens were preheated to 436 K 
(32S(> F) . The cold mold was placed in the press and contact 
pressure was applied and maintained for 3 additional minutes after 
the thermocouple indicated 311 K (100<> f) . Then the pressure 
was increased gradually over a period of 2.5 minutes to 0.34 MPa 
(50 psi) . Press temperature and pressure were then held constant 
for an hour at which time the power was turned off and the 
laminate allowed to cool under pressure. The laminate was 
postcured in the mold under a weight of 4.54 kg (10 lb) for 16 
hours at 411 K (280° F) . 

Specimen Preparation, Instrumentation, and Testing 

The laminates were cut into test specimens 2.54 cm (.'-in.) 
wide using a dieunond cutting wheel. Specimen length was .0.5 cm 
(12-;n.). Specimen ends were reinforced with 5.1 cm (2-in.) long 
adhesively bonded fiberglass end tabs leaving gage lengths of 20 
cm (8-in.). Each tensile specimen was instrumented with one 
120-ohm, 60° delta rosette strain gage located at midlength for 
specimens that were to be tested monotonically to fracture. So as 
not to ir'terfere with subsequent frequency testing, the gages were 
located at about 1.2 cm (0.5 in.) from one end of the test section 
on the specimens that were to have damage induced via cyclic load- 
ing. 

A specimen of each material was loaded to fracture in a hy- 
draulically actuated universal testing machine. Testing was in- 
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crcMiital to facilitate periodic recording of strain gage data. 
Stress-strain curves were also plotted concurrently with the test- 
ing to observe changes in slope of these curves that eight indi- 
cate daeage. Based on inforeation obtained froe the destructive 
testing of these two specieens. other specieens were subjected to 
cyclic loading at stresses large enough to induce daeage but not 
sufficiently large to cause fracture. Daeage thus achieveo was 
revealed by ultrasonic C-scans of specieens treated as described 
above. Specieens 10.2 ce (4-in.) long were taken froe selected 
areas of these tested tensile specieens and froe untested eater ial 
for detereination of vibration frequencies ana deeping factors via 
low velocity iepact excitation tests. 

Measureeent of Frequencies and Daeping Factors 

The vibration frequencies and the daeping factors were deter- 
eined by a special procedure which was developed at the Oniversity 
of Cincinnati under contract to NASA Lewis Research Center. The 
tests and data reduction were conducted at the University of 
Cincinnati as a part of this saee contract, in this procedure the 
test specisien is considered to be a free beaa. Briefly, then, the 
procedure is as follows. The bean, suspended vertically by n«eans 
of low swss silk threads bonded to the heasi eoge, was excited with 
an isgwct froe a saall ball bearing and the response was Mea- 
sured with a wide- frequency response aicrophone. n>ere were no 
transducers in contact with the bean. There was no swasurable 
loading of the be as. A zoom transfora was used with a LaPlace 
transfoca algor ithn (Fast Fourier Transfora) to compute both the 
natural frequency and daaping factor. The beaa was tested with a 
nuaber of identical isguicts at different tines. The aeasurenents 
were reproducible within a naxinua standard deviation of 0.0013 
percent for natural frequencies and 1.3 percent for aodulus of 
elasticity. This error, of course, can be reduced by averaging a 
nuaber of impacts by a 1/ ^n factor. 

The zoom transform was used to achieve a small delta frequency 
and the Laplace transform was used to compute the natural fre- 
quencies from the zoomed data. 

Since the beam was small, it was not practical to mount a 
transducer on the bean. Therefore, the mode shape information was 
determined by mounting a small known mass on the beam ana measur- 
ing the change in natural frequency due to the mass. The change 
in the natural frequency can be used to compute the effective mass 
of the beam's modes of vibration at the point where the mass has 
been added. By moving the mass over the surface of tne beam, the 
mode shape can be measured. 

The effective mass of the mode of vibration can be determined 
by the following relationship: 



where Mg » 

effective mass 

Ma = 

added mass 

Wi 

initial natural frequency 


natural frequency with added mass 


THEORETICAL PART 


The theoretical part of the investigation consisted of 
laminate analysis and finite element analysis via NASTRAN. 
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l.uiinate Analysis 


Laainate analysis was usea to predict the laainate aoduli, 
ply residual stresses, ply stresses at fracture and the aate- 
rial properties needed for input to NASTRAN, The actual calcu*' 
latlons were perforaed usin 9 the coaposite aecnanics coaputer 
code^. The daaage induced in the laainate by the residual 
stresses, the preloading and the cyclic loading was siaulated 
indirectly using an equivalent void voluae ratio in the coa- 
puter code. In essence this approach considers the laainate to 
have unifora daaage distribution which is acceptable for fre- 
quency calculations but iiut for local stresses. In addition 
the effects of possible shear coupling were siaulateo by per- 
aitting the laainate to be slightly anisotropic in its plane. 

Finite Eleaent Analysis 

The finite eleaent analyses were perforaed to predict the 
free-beaa vibration frequencies of the speciaens ano the ef- 
fects of low level daaage, due to factors aentioned previously, 
on these frequencies. Additional analyses were perforaed to 
assess the effects on the frequencies of thruugh-the-thickness 
defects in the aiddle of the speciaen. The finite eleaent 
aodel of the specimen that was used in NASTRAN^ is snown in 
Figure I. It consists of 297 nodes, 2S6 quadrilateral aniso- 
tropic plate (CQUA02) eleaents, and 14F3 degrees of freedoa 
(OOF) . This type of aodel was considered reasonable to siau- 
late the effects of low level damage on the free ’ ibration re- 
sponse of these specimens. The through-the-thickuess defects 
were simulated by removing two elements (124 ano 125) from the 
center of toe specimen, 

RESULTS, COMPARISONS ANO DISCUSSION 

The experimental and predicted results obtained in this 
investigation for the various properties described previously are 
summarized, compared and discussed in this section. Comparisons 
between experimental and predicted results are also summarized. 

Stress-Strain Curves 

Typical monotonic ano cyclic stresr--strain curves of the 
specimens tested are shown in Figure 2. These curves indicate 
changes in the initial and tangent slopes tor botn the E-G/E and 
the T-75/E composites under both monotonic and cyclic loading. 
This indicates thac both loadings nave inouced damage in the 
specimens. Note also the "apparent hysteresis" in the cyclic 
loading curves. The corresponding I'oisson's strain curves ace 
shown in Figure J. As can be se ;n in this figure, the Poisson's 
strain curves show similar behavior as rh“ stress-strain curves, 
C-scan traces for tho specimens are shown in Figure 4. These 
traces- show conside’.ubly more damage for the T-73/E composite 
specim'ins than for che E-G/E specimens. Photographs of typical 
specimens tested including some of those tested to fracture, ace 
showr in Figure 1. 


Ply Stresses 

Unidirectional and angleplied composite properties of the 
specimens required for the various analyses were predicted using 
the computer code^. These properties ace summarized in table 
I. Limited available experimental data for these specimens ac!; 
also shown for comparison. 
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As can be seen the predicted values are in qood aqreeaent with 
the available expertaental data (within 10 percent for aoduli). 

The calculated transverse resiaual stresses in a given ply of 
the angleplied lasiinates were 25 MPa (1^6 iisi) for the E-G/B and 
33 NPa (4.8 tsi) for the T-7S/E. The corresponding longitudinal 
stresses have the sane nagnitude but opposite sign, while the in- 
tralaeinar shear stresses are ‘zero". These types ot resiaual 
stress states are characteristic in quasi-isotropic angleplied 
IsBinates. 

Calculated ply stresses at fracture loads, obtained using the 
properties shown in table 1 and including residual stresses, are 
sussnarized in table II. The combined-stress fcilure-criterion 
values for the plies are given in the column *Hargin-ot-Safety* 
(Most. Negative values in this column indicate ply failure. Typ- 
ical uniaxial ply strengths (fracture stresses) are also given in 
this table. Based on the negative HOS values, the * plies and the 
90° plies failed in both the £-G/E and the T-75/E laminates. 
Coaqiaring ply stresses with the corresponding ply strengths, it is 
seen that the M5° and 90° plies failed in transverse ten- 
sion. About one-fourth of the transverse ply stress in the C G/C 
laminate is due to residual stress as is about 60 percent in the 
T-75/E laminate. It is worth noting that the ply transverse re- 
sidual stresses in the T-75/E laminates 33 MPa ((4.8 ksi)) are 
greater than the corresponding ply strength (29 MPa (4.0 ksi)) and 
cause transply cracks in the laminate as shown in Figure b. The 
transverse ply stresses, at the cyclic stress shown in the ultra- 
sonic C-scan ifig. 4c), are respectively: 85 MPa (12.3 ksi) in the 
♦45° plies of the E-G/E laminate and 41 MPa (6.0 ksi) in the 
+45° plies of the T-75/E. Both of these are greater than the 
cor resi>onding transverse ply strengths in Table Ilb and, there- 
fore, induce damage (transply c.acks) in these laminates. The 
corresponding transverse ply stre:.ses in the 90° plies are about 
twice as large and the damage will be more severe in these plies. 
Recall that the lAagnitude of the cy :lic stress was selecteo based 
on the changes in the ■slope of the £ tress-strain curves shown in 
Figure 2. The calculated ply stresses just described show which 
ply stresses induce the damage. 

The above discussion demonstrate directly and indirectly that 
the cyclic load used in the invest) ■; tion induced damage in the 
0°, +45°, and 90° plies in both E-G,- and T-75/E laminates. 

Vibration Modes and Damping Factors 

The measured results, obtained from the special procedure de- 
scribed previously, for free vibration frequencies and damping 
factors were determined from frequency response spectrums such as 
those shown in Figure 7. The reduced data are summarized in Table 
HI. Note the results shown are normalized with respect to those 
for the undamaged specimens for ease of comparison. 

It can be seen in Table III that the frequency decreases with 
increased load (damage) for the E-G/E laminates as was expected. 
The decreases in frequency range from about 3 percent to 6 per- 
cent. The damping factor increases ranging from about 1 percent 
to about 29 percent for the lower load case and from about 3 per- 
cent to 18 percent for the higher load case. A., exception to this 
is the damping factor for the second frequency which decreases for 
the lower load case and increases for the higher load case. This 
indicates some kind of an anomaly that may be attributed, in part, 
to the degree of coupling between bending and torsion. The con- 
clusions to be drawn from these results are; 

1. The frequency of E-G/E laminates decreases in the presence of 
damage although this decrease may be small and may be shadowed oy 
the effects of other factors; and (2) the damping factor for E-G/E 
laminates increases in general with no dear correlation between 
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damping and damage. 

The results for the T-7S/E laminates in Tatie 111 show pro- 
gressive rncreases in frequencies with ihcreased load (damage) and 
a mixture of increases and decreases for damping. These mry be 
caused, in part, by the high stiffness which minimizes the vibra- 
tion amplituder. The conclusion from these results is that fre- 
quencies or damping factor may not be sensitive discriminators for 
damage assessment ir T-7S/E laminates. 

The experimentally determined frequencies and those predicted 
using NASTRAN are summarized in table IV for comparison. The fol- 
lowing points are worthy of note. For the E-G/E laminate: 

1. Two more frequencies were predicted than were determinea by 
measurement; 

2. The predicted frequencies are in good agreement with the cor- 
responding measured ones (within S-percent) for the undamaged 
case (compare column 1 with column 5); 

3. Twenty percent void equivalent in all the elements was needed 
to simulate the damage induced by the loads (column 8) ; 

4. Simulation using shear coupling (column 6) and damage in se- 
lected regions in the specimen (column 7) , based on the 
C-scans (fig. 4c), is not sensitive enough to predict the 
measured reiuction in frequencies (column 2 or 3); 

5. The tree frequencies are insensitive to local danage 
(through-the-thickness defect at the specimen center, with 
defect- length/specimen width equal to 1/4) , (column 9) . 

For the T-75/E laminate; 

1. Three frequencies were determined by measurement which were 
not predicted using NASTRAN and two frequencies were predicted 
which were not determined by measurement 

2. The predicted frequencies are within abouL 3 percent of the 
corresponding measured ones for the undamaged specimens (com- 
pare columns 1 with 5); 

3. The measured frequencies are not sensitive to load ard resid- 
ual stress damage as was already mentioned in the previous 
section; 

4. Shear coupling and partial damage have negligible effects on 
the predicted frequencies (less than 2 percent, columns and 
7); 

5. A uniform damage simulated by 20 percent voids in all the’ ele- 
ments reduces tne predicted frequencies by about 11 percent 
(column C) ; 

6. Concentrated local damage has negligible effect on the free 
vibration frequencies (column 9) . 

It is apparent from the above discussion that both experi- 
mental and analytical methods may be needed to determine the vari- 
ous modes in a vibration analysis. It is also apparent that the 
experimental procedure described herein, low-velocity impact 
coupled with the Fast-Four ier Transform and packaged in a minicom- 
puter, can be a convenient procedure for assessing low-level dam- 
age on the dynamic response of fiber composite angleplied 
laminates. 

The first five vibration mode shapes corresponding to the vi- 
bration frequencies column 9, Table IV are shown in Figure 8. As 
can be seen in this figure the mode shapes are primarily bending 
and torsion. It is noted that frequencies for the first, second 
and third bending modes are not integer multiples of each other. 
Neither are the two torsional frequencies. This means that each 
of the mode shapes in Figure 3 has some degree of coupling between 
bending and torsion. Typical mode shapes plotted from the data 
obtained from the experimental procedure are shown in Figure 9 for 
the first and second bending modes. As can be seen by inspection 
thece modes are in excellent qualitative agreement with the cor- 
responding predicted ones in Figure 6. 
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CONCL.Ui>IONS 


The major conciusxons ot an investigation to assess the sen- 
sitivity of the dyniunic response (free vioration frequencies ana 
damping factors) of fiber composite angleplied laminates to low- 
level damage induced by load (monotonic and cyclic) and resioual 
stress are as fellowis: 

1. Monotonic load, cyclic load and/or residual stress induced 
damage may reduce the free vibration frequencies and may in- 
crease the damping factor in low modulus fiber comp^isites such 
as E-Glass/Epoxy (E-G/E) angleplied laminates, while it may 
have negligiule effect, if any, in ultra high modulus fiber 
composites such aS Thornel-75/Epoxy (T-75/E) angleplied lamin- 
ates. 

2. Proper simulation of load induced damage for analysis may re 
quire considerable judgment and intuition. 

3. Free vibration tests may not he sufficiently sensitive to as- 
sess concentrated local damage in angleplied composite lamin- 
ates. 

4. Both experimental and analytical methods may be needed to oe- 
termine the various modes in vibration analyses. 

5. Dynamic response determined from low-veiocity impact coupled 
with the Fast-Four ier-Transform and packageo in a minicomputer 
can be a convenient proceoure for assessing low-level damage 
in fiber composite angleplied laminates. 
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^Fitted slope. 

CoBvcrslon fsetors: 10* psl - O.SMt GPs 
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TA?LE n. - CALCULATED PLY STRESSES AT FRACTl RE. MCLUDING RESIDC^ STRESSES, 
ANDCOHKESPONUNG UNIAXIAL PLY STRENGTHS 
(a) Ply jtressos 


Composite 

Frai'ture 

Ply stresses, ksi 


Bt-ength, 


145® 

Plies 



90° 

Plies 


0® Piles 


ksi 

Longl- 

tudtrsl 

Trans- 

verse 

9iesr 

SlargiB 
of safety 

U>ngi- 

tudinai 

Trans- 

verae 

Shear 

'tunu 

of safety 

Longl- 

btdiAal 

Trans- 

verse 

jbear 

®Msrglii 
of safe^ 

E-glMs/epoxy 

59.6 

41.3 

16.3 

‘16.9 

-38 

-IT. 9 

31.7 

0 

-63 

100 

5.0 

0 

-0.36 

T-TS/epox!' 

44.8 

38.6 

C.O 

1.6 

-.S3 

-25.T 

7.0 

0 

•2.7 

103 

4.9 

0 

.67 


t'nlajUal ply strengths 


Composite 

1 Uniaxial ply strength, ksi | 


LongltudlAal tensile 

Longitudinal compressive 

Transverse tensile 

Transverse compressive 

In-plane shear 

E-glass/epoxy 

170 

97 

4 

17 

5 

T-75/epoxy 

142 

126 

4 

19.8 

7.1 


‘Combined sircss ftlfure. 

Cooversioo factor: CPs • ksIv 0. 66»47xl0' 






TABLE tv. - FREE VIBRATIDN FREQUENCIES OF UND/>MAOED ANP DAMAGFD COMPOFITE8 


Composite 

Fre- 

quency 

oraer 

Frequency, Hz 


EjqMrlmental 

I Analytical 

1 

2 1 3 

4 

5 

6 

7 



E^glaaa/epoxy 


Undamaged 

After loading to 


Undamaged, 
no ahear 
coupling 

Undamaged, 
with shear 
coupling 

PartliJ 
damage with 
shear coupling* 

All elements 
j damaged^ 

Slotted^ 

0.193 GPa 
<28.2 kel) 

0.333 GPa 
(48.3 ksi) 


1 

630.1 

2246.1 

781.2 

3197.3 

781.2 

2246.1 


815.3 

2264.3 

2432.4 

4469.4 
4990.3 
7426.8 

813.8 

22.59.8 

2427.5 

4460.0 

4980.1 
7 .10.3 

803.9 

2231.6 

2399.7 
4:^6. 6 
4923. S 
7374.5 

738.2 

2051.0 

2187.6 

4049.7 

4490.0 
6730.4 

802.4 

2259.4 

2404.5 

4427.4 

4982.4 

4296.9 

4199.2 

4150.0 

7080 

6933.6 

6884.8 

T-75/epoxy 

■ 

Undamaged 


After loading to 

Undamaged, 
with shear 
coupling 

Partial 
dan age with 
shear coupling 

All elements 
, damaged 

Slotted 

0.036 GPa 
(S.2 kzi) 

0.299 GPa 
(43.4 kal) 

0.309 GPa 
(44.8 ksM 

nc shear 
coup'<ng 

1 

1035.1 

1051.7 

1060.9 

llOf.O 






2877.9 

4082.6 

5643.7 
80P9.6 
9254.5 

2912.4 

4092.4 
5692.3 
8124.2 
9316.8 

2962.9 

4046.2 
5730. S 

8000.9 

9317.2 

3018.6 
4177.1 
6002.4 
8283.9 

9814.6 

1395.9 

1385.9 

3907.0 

5453.9 
7819.4 

1366.4 

1241.2 

1368.9 

3906.6 

5399.3 

7763.3 

3940.2 

5490,6 

7887.4 

3852.8 

5385.7 

7661.4 

3501.8 

4879.8 
7008.6 

13550 

13658 

13660 

14295 

11181.6 

13252.5 

11099.6 

13148.9 

10962.0 

11064.0 

5933,4 

11776.4 

11104.7 

13134.8 


Wli input for element. 57-64, 121-128, m,d 195-200 b end on 0.2 void volume frKtlon. 
^oAill inpu* for all elementa baa«<i on 0.2 void volume fracUon. 
fflot rq>reeented by removal of elomenta 124 and 125. 

















































































4CYCUSBEP«0ll5AND48k$i 3 CYCLES BEIWEIN 5 AND O ksi 

to CYCLED BaOW FRACTURE POINT. 

Figure A - C - Scans (4 untested composites and o( tested specimens H ksi MPa). 
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!b) T-^^/EPOXY. 
figure 5. - Tested spetimens. 
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Rqure 6. - Photomicrogrph shewing transply cracks in [+45/90o/0,], Thornel 
75^epoxy angleplicd laminate. c ca 
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